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Abstract. We show the generic existence of metastable massive gravitons in the four- 
dimensional core of self-gravitating hypermonopoles in any number of infinite-volume 
extra-dimensions. Confinement is observed for Higgs and gauge bosons couplings of the 
order unity. Provided these resonances are light enough, they may realise the Dvali- 
Gabadadze-Porrati mechanism by inducing a four-dimensional gravity law on some 
intermediate length scales. The effective four-dimensional Planck mass is shown to be 
proportional to a negative power of the graviton mass. As a result, requiring gravity 
to be four-dimensional on cosmological length scales may solve the mass hierarchy 
problem. 
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1. Introduction 

The old and more recent interest in the existence of space-hke extra-dimensions has led 
to three main ways of accommodating their presence with an apparent four- dimensional 
world [HEIIB]. In String Theory, extra-dimensions are compactified and can be probed 
only at an energy scale exceeding their inverse radius. The non-observation of the 
associated Kaluza-Klein (KK) excitations at colliders have pushed such a scale above 
the TeV. However, motivated by orbifold constructions, it has been realised that the 
typical size of the extra-dimensions probed by gravity could be much larger than the one 
felt by the gauge (and matter) fields [HEIE]. In this picture, our universe could be a four- 
dimensional brane on which the Standard Model particles are confined, embedded in a 
higher- dimensional bulk [3 [U [9l [TOl [11] . Still, gravity is tested to be four- dimensional 
on length-scales ranging from the micrometers to the cosmological distances [III [13]. 
A micrometer is nevertheless much larger than a TeV~^ while the current cosmic 
acceleration might be the signature that gravity is actually no longer standard on the 
largest length scales [H]. These motivations have led to two other gravity confinement 
mechanisms. In the Randall-Sundrum (RS) models, the extra-dimensions felt by gravity 
are non-compact but still of finite volume such that deviations from four- dimensional 
physics appear only at lengths smaller than their curvature radius [15]. In the Dvali- 
Gabadadze-Porrati (DGP) models, the extra-dimensions are non-compact and of infinite 
volume. Gravity can be observed as four-dimensional provided the graviton flux between 
two masses is prevented to leak into the extra-dimensions. In the original DGP 
model, this is obtained by adding a (quantum induced) four-dimensional Einstein- 
Hilbert term on the brane [161 HZ] ■ Although it has been argued that some instabilities 
should show up in the original approach [18], [I9] , stable extensions have been proposed 
since [2ni[2ll[22l[23l[2ll[25l[26l[27]. The so-called regularised DGP models assume 
an ad-hoc varying gravitational coupling constant along the extra-dimensions such that 
gravitons are reflected back onto the brane in the same way as a varying dielectric 
constant reflects photons. Assuming the gravity action to be 



where stands for the number of extra-dimension, R is the Ricci scalar and g is 
the determinant of the metric tensor, it has been shown that 4D-like gravity can be 
recovered provided the function J^{X) is peaked enough in a narrow region around the 
brane [281 EH] . However, the existence of a similar mechanism in a well-defined physical 
framework is a non-trivial issue. Indeed, the role of the function J-" in Eq. ([1]) could 
be played by g{X) if there are non- vanishing stress tensor sources in the bulk |30j . 
A varying Planck mass can also be obtained by considering a dilaton ip that would 
condense on the brane such that J-'{X) = exp[ijj{X)] has the required profile. In any 
case, none of these fields would be independent and their equations of motion have to 
be solved in a given field theoretical setup to address this question. 




(1) 
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This has been done in Refs. [311 [22] in a scalar-tensor theory of gravity by 
considering our brane to be a self-gravitating topological defect [33l [Ml ISSl [361 [37] . 
The core is a flat four-dimensional spacetime supposed to be our universe while the 
non-trivial defect-forming field configurations curve the extra-dimensions. Such a defect 
can be formed by the spontaneous breakdown of a local symmetry in the bulk whose 
stress-tensor ends up being non-vanishing only in a region defining the brane thickness. 
Typically, it is given by the Compton wavelength of the various defect-forming fields. 
The extra-dimensional spacetime is asymptotically flat and of infinite volume, as in 
the DGP model. In six dimensions, this can be obtained by breaking a local U{1) 
symmetry to form a hyperstring whose core has three spatial dimensions (there is rig = 2 
extra-dimensions). The seven-dimensional version is a 't Hooft-Polyakov hypermonopole 
obtained by breaking an S0{3) symmetry in ne = 3 extra-dimensions. As shown in 
Ref. |3]J|, reflecting gravitons onto the defect core requires a violation of the positivity 
energy conditions in General Relativity [38]. This happens to be impossible for the 
hyperstring unless one adds a source of negative energy in the bulk. As a matter of fact, 
a negative cosmological constant does the job and the model becomes of the RS type, 
with finite volume extra-dimensions and a severe fine-tuning problem [391 [10]. It is well 
known that a constant positive curvature term in the Einstein equations behaves like a 
perfect fluid with a negative equation of state parameter and can therefore mimic matter 
with negative pressure [J^. Ref. [32] has shown that this mechanism is indeed at work 
inside the seven- dimensional hypermonopole: assuming isotropy, the positive curvature 
of the two-dimensional orthoradial extra-dimensions acts as a potential barrier for the 
propagation of gravitons. These become resonant and massive. The next question is 
therefore to determine if curving more than two orthoradial extra-dimensions still allows 
for a similar gravity confinement in higher dimensional spacetimes. 

The present article is devoted to this issue. In the following, we show that the DGP- 
like graviton confinement mechanism by curvature effects is indeed a generic feature 
of self-gravitating hypermonopoles and occurs in any asymptotically flat spacetime of 
strictly more than six dimensions. More precisely, in 4 + ne dimensions, hypermonopoles 
can be formed by the spontaneous breakdown of an 5*0 (rig) symmetry to SO{ne — 1). 
Moreover, in order to allow for a varying Planck mass in the bulk, gravity is supposed to 
be of the scalar-tensor type, ip being the dilaton. After having derived and numerically 
solved the equations of motion in Sec. [21 we show that the extra-dimensions are of 
infinite volume and asymptotically flat while being strongly positively curved in an 
intermediate region. Such a conflguration is obtained without flne-tuning and naturally 
occurs for values of the fleld coupling constant of the order unity. In Sec. [HI we solve 
for the propagation of spin-two fluctuations along the brane and flnd strongly peaked 
massive metastable resonances. We then illustrate how these resonances may realise 
the DGP mechanism by deriving the resulting Newtonian potential on the brane: it 
remains (4 -|- rig) -dimensional at small and large distances but behaves as dimensional 
in an intermediate range with d < 4 + ne. Finally, we discuss the mass hierarchy 
problem and show that the four-dimensional effective Newton constant is proportional 
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to a positive power of the mass nig of the hghtest associated graviton resonances. This 
property is analogous to the existence of a cross-over distance at which gravity becomes 
iV-dimensional in the DGP regularised models [281 EH]. In terms of the Planck mass, 
we find in Eq. (175|) that 

^Pi ^ ' (2) 

where M stands for the A^- dimensional Planck mass. Since four-dimensional gravity 
requires extremely light resonances, our mechanism necessarily implies a very small four- 
dimensional gravity coupling constant. Finally, concerning the smallest length scales, 
we show that the distance under which gravity is again A^- dimensional depends on the 
gravitational redshift induced by the hypermonopole forming fields. It can be made 
arbitrarily small, independently of the graviton masses, provided the Higgs and gauge 
bosons have masses close to the Planck mass M in dimensions. 



2. Hypermonopoles of any dimension 

In this section, we assume our universe to be the four-dimensional core of a 
hypermonopole living in = 4 + Ue dimensions. This topological defect can be 
formed by spontaneously breaking an SO{nc) symmetry to SO{nc — !)• We impose 
this symmetry to be local such that the defect does not exhibit long range interactions 
and has a localised stress tensor allowing asymptotically flat extra-dimensions. As 
mentioned in the introduction, gravity is assumed to be of the scalar-tensor type to 
allow for a varying effective Planck mass along the extra-dimensions. In the Jordan 
frame, the action describing this system reads 



X 



+ V\9 



Z o 



x, (3) 



where the SO{ne) Higgs field is an rig- dimensional vector $ = {0''} with a G 
{1, 2, . . . , ne}. Gauge invariance under local SO^ric) transformations is ensured by 
defining the covariant derivatives as 

Va^ = Oa^ - qCa^, (4) 

where q is the Higgs charge and Ca = {{Ca}"'''} are the gauge field matrices. The 
associated field strength tensor matrices are given by 

Uab = OaCb - ObCa -q[CA, Cb] . (5) 
The last term in Eq. ([3]) is the Higgs potential V{^) and breaks SO{nc) to SOijic — 1) 
such that the topology of the vacuum manifold is the same as the n-sphere 5" = 
SOijif.) / SOijic — 1), where we have defined 

n = rie — 1. (6) 
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Since the n*"^ homotopy group '7r„(S'") = Z is non-trivial, we expect the formation 
of hypermonopoles mapping SO{nf.) to the extra-dimensions [12]. In Eq. ([3]), the 
quantity U{ip) encodes the dilaton potential in the Jordan frame. For simplicity, we 
choose the dilaton to be a free massive particle of mass in the Einstein framqj such 
that 

Uiij) = m2e2^/("+=^) . (7) 

One can check that the above equations reduce to the hyperstring of Ref. [31] for = 2 
and to the 't Hooft-Polyakov hypermonopole of Ref. [32] when Ue = 3. 

In the following, after having introduced our Ansatz for the field profiles, we derive 
and solve the equations of motion assuming isotropic extra-dimensions. 



2.1. Equations of motion 

2.1.1. Gravity sector Varying the action with respect to the metric tensor and the 
dilaton gives the Einstein- Jordan equations 

Gab = e-^^'TAB + ^Oa^Ob^ - \gAB [39x^9^^ 

+U{ilj)] + WaOb^ - gAB^iP, (8) 

where □ = V xd'^ ■, Gab is the A^-dimensional Einstein tensor and Tab the matter stress 
tensor 

Tab = Va^ ■ - ^g^'^MnAc'HBD) + QabC^, (10) 
where C^a is the Higgs and gauge field Lagrangian. 



2.1.2. Matter sector The variations of Eq. ([3]) with respect to the Higgs and gauge 
fields gives the Klein-Gordon and Maxwell-like equations: 



and 



dV 

- qg^"" [(V^Cb$) + CbVa^] - = 



(9$ 



{VAn^^'f +qg^'' iVA^r<i^' - {Va^ 



AB\<^b 



+ q {n^"" -Ca-Ca- n^"") 



0. 



(11) 

(12) 
(13) 



2.1.3. Metric and field Ansatz Respecting the hyperspherical static symmetry in the 
rZe extra-dimensions and Poincare invariance along the brane gives the metric 

ds2 = ^^^dx^dx^ = e'^'^^r/^^dx'^dx'^ + dr^ + w^(r)d(]^_ ^^4) 



X in which the scalar and tensor degrees of freedom are decoupled 
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n 

dQl = Y,'^.,{e,^,)del (15) 

1=1 

with Ti = 1, T2 = sin^(6'i) and 

i-l 

TiiOj^i) = T2{9i) II sm\9,) for i = 3,...,n. (16) 
i=2 

For a defect configuration, the Higgs field vanishes in the core whereas it asymptotically 
recovers its vacuum expectation value. Enforcing the spacetime symmetries, we assume 
a radial field such that 

r = ^ir)-, (17) 
r 

with ip{0) = and ip{r) — > v when r — )■ 00. In Eq. ( IT7|) . = 6abX°'x^ where the {x""} 
stands for the Cartesian coordinates defined by 

=rcos9i, = r sin 6*1 cos 6*2, (18) 
x"' = r sin 61 sin 62 . . . sin 6'„_i cos On, 

x"" = r sin 61 sin 62 ■ ■ ■ sin ^„_i sin 6'„. (19) 

We also assume that the dilaton depends only on the radial coordinate 

ip = ip{r) . (20) 

Our Ansatz for the gauge field is the generalisation of the 't Hooft-Polyakov 
configuration [lH HI], with a unity winding number. Requiring the stress tensor to 
vanish at infinity imposes the covariant derivative to vanish. From Eqs. (jl]) and ( I17p . 

one gets 

(C. r'' - ^ - ^^'^ f^x" - ^xA (21) 

where the dimensionless function Q{0) = 1 for regularity in the core and Q{r) — at 
infinity. All the other Ca are vanishing. Under such an Ansatz, observing that 



Tr{CoA^) = -2 



l-Q(r) 



2 



, , 9e.e,, (22) 
quj[r) 

the gravity and matter sector equations considerably simplify and we write down only 
the final result in the next section (see the appendix for some intermediate steps). 



2.1.4. Dimensionless equations For convenience, we introduce the following 
dimensionless quantities. The radial distance can be expressed in unit of the Higgs 
Compton wavelength such that 

p = niiir = VXv r, (23) 
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where is the mass of the Higgs boson. Similarly, the dimensionless angular metric 
coefficient and Higgs field are defined by 

zu = miyU, f = — . (24) 

The gravity, Higgs and gauge coupling constants account for three dimensionless 
parameters in the equations of motion ([H]) to fll2p which can be recast as 

..^4 = 4. = (25) 



mb being the mass of the gauge bosons. After some rather long algebra, a dot denoting 
differentiation with respect to p, the dimensionless equations of motion in the gravity 
sector read 

3 3 2 ^ / ~ 1 3 w 

= -ae^£:-V^-^^^-0a + n|)^-|, (26) 

3.2 / X ri?^ — 1 „ . 
-a +n{n-l) ^ +2na— 

= ae-^V + - (2& + n-] (27) 
2 \ zu / Z 

5 r. , .zu , . , . ijp' — 1 , . zu 

2a + -d^ + {n- 1)— + (n - l){n - 2) ^ + 2{n - l)a— 

2 zxj 2zu zu 



, , zu 
2a+{n- 1) — 
zu 



i'-j, (28) 

^ +-^^+^2* + n^j^ = -^t/ + ^-flj. (29) 

The dimensionless dilaton potential is f/ = U / m\ and the dimensionless Ricci scalar 
R = R/m\ stands for 

• 2 1 

R = -Aa -hb^ -2n 72(72 - 1) And— . (30) 

zu uu zu 

The quantities £ and V are respectively the energy density and pressure generated by 
the Higgs and gauge fields along the radial extra-dimension T 



z 2ezu 2zu 4:ezu 



^ = V + + w^fQ" + '^V._/' (i - Q')' + (31) 

P = V + S - ^f'Q' - - ^')' - ^' (32) 

z Zezu Zzu 4:ezu 

where V{f) stands for the dimensionless Higgs potential 

V{f)^l{l-f)\ (33) 
§ We also have £ — -Cm/im'^v'^). 
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The quantity £± appearing in Eq. (1251) is the energy density along the orthoradial 
directions and reads 



^± = V + o 2 + ^T^/^ + ^ 7^^ ^(1 - Q ) + y- (34) 

The dynamical equations for / and Q stem from Eqs. flTT]) and flT2|) 

/-n^- - — 

/ zu dj 



f+ {2a + n^]f-n^-^ = 0, (35) 



2cr + {n- 2) — 

Ct7 



Tl — 1 

Q + — ^(g-g')-e/'g = o. (36) 



These equations match with those of the six- dimensional hyperstring and seven- 
dimensional hypermonopole derived in Refs. [211 E2] • Let us however notice the presence 
of new terms for n > 2 in the orthoradial equation (!28|) as well as in the associated stress 
energy in Eq. fl34p . It is also worth remarking from Eq. fl34p that for n = 2 or 3, the 
gauge field is generating a negative potential in the orthoradial extra-dimensions, which 
becomes again positive for n > 4. 

2.2. Background fields and geometry 

2.2.1. Boundary conditions The boundary conditions for the metric coefficients 
and fields are fixed by requiring regularity in the core and a Dirac hypermonopole 
configuration asymptotically. As already discussed, we look for asymptotically fiat 
spacetime, i.e. 

lim /(p) = 1, hm Q{p) = 0, 

lim a{p) = 0, lim = 1, lim ^(p) = 0. (37) 

p— >oo p— >-oo p p— >oo 

Notice that a could be shifted by a constant value since all the equations of motion 
depend only on &: this refiects the expected invariance with respect to a rescaling of 
the internal brane coordinates x'^. We have also chosen the dilaton to vanish at infinity 
since this minimises its potential energy and t/) = is an exact solution of Eq. (129!) for 
^ = 0. One can check that this last condition is indeed asymptotically fulfilled with the 
limits of Eq. ( !37l) . Let us notice that the metric far from the core is not a generalisation 
of the conical fiat metric existing around a cosmic string. As can be checked in Eq. f l5U]l . 
as soon as n{n — 1) 7^ one has w = p, i.e. there is no missing angle (for & = 0). 

In the hypermonopole core, the SO{ne) symmetry should be restored and the 
spacetime geometry has to be regular. As a result, the fields satisfy 

lim/(p)=0, limQ(p) = l, 

p— >-0 p— >-0 

lim (t(p) = 0, lim^^ = 1, lim^(p) = 0. (38) 

p->0 p—^O p p->0 
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Figure 1. Field and metric profiles forming an eight-dimensional (n = 3) 
hypermonopole for a — 2.60, e — 0.25 and /3 = 1.00 (top). The space-time is fiat 
in the core and at infinity but strongly curved in the intermediate region. The dilaton 
condenses in the core and traces the energy density distribution. See also Fig. [2j 

2.2.2. Solutions From Eqs. fl371) and fl38l) . we have ten boundary conditions to solve 
the ten-dimensional first order non-linear differential system that can be obtained from 
the second order Eqs. ([26]), ([28]), ([29]), ([35]) and ([36]). Notice that Eq. ([27]) is not 
included since this is a constraint equation and is redundant with the previous set, 
up to a constant which is fixed once the boundary conditions are specified. Finding 
numerical solutions of this system is non-trivial and we have used the conditioning mesh 
methods implemented in Ref. [15]. We have first checked our numerical implementation 
by recovering the n = 1 and n = 2 solutions of Refs. [311 [32] before solving the system 
for n > 2. Despite the new terms appearing in the equations of motion, we have 
found hypermonopole solutions for any tested value of n. All of them exhibit similar 
patterns than those found in seven- dimensions. For coupling constant of order unity, 
the spacetime is strongly curved in an intermediate region where the field derivatives 
are non-vanishing, and in particular zu remains almost stationary with respect to p. 
As a result, the hypersurface of the n-sphere of radius p becomes constant and the 
extra-dimensions are cylindrically shaped. At further distances, w ^ p again and the 
spacetime becomes flat. 

In Fig. [H we have represented the field profiles obtained for n = 3, i.e. in eight 
spacetime dimensions. The dilaton condenses in the core as a scalar gravity field 
which passively follows the stress energy distribution. The metric factor a traces the 
gravitational redshift between the core and the asymptotic spacetime and has been 
represented with uir) in Fig. [21 Finally, the Higgs and gauge fields are typical of a 
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5 10 15 20 25 30 

Figure 2. 5" hypersurface in tlic strongly curved region: uj{r) remains almost constant 
and the extra-dimensions becomes cylindrically shaped. The gravitational redshift in 
the core is given by e'^'-"-'. The parameters are the same as in Fig. [TJ 

topological defect configuration. By choosing e = 0.25, the gauge bosons are twice 
lighter than the Higgs boson and condense within a larger extra-dimensional radius. 
The shift in the condensation radius of the Higgs and gauge field produces the step 
observed in the dilaton profile (see Fig. [1]). 

2.2.3. Dependence in the number of dimensions As discussed in the beginning of this 
section, the hypermonopole-forming fields exhibit a similar behaviour for all values of n. 
In fact, this can be understood from the equations of motion f l26|) to fl36|) . The number 
of dimensions enters these equations through n at two levels. 

First, it changes the coupling between the metric coefficients w and a, as well 
as w and the other fields. For a given xxj{r) profile, one expects all the fields to 
be more sensitive to it when n becomes large. Let us emphasise that this effect is 
purely geometrical since it persists when a — ?► 0. In Fig. |3l we have represented the 
hypermonopole solutions obtained for an assumed generic set of parameters a = 0.1, 
e = 0.2, /3 = 1 and in various dimensions ranging from seven (n = 2) to ten (n = 5). 
The upper panel of this figure illustrates the above-mentioned effect. For larger values 
of n, the gravitational redshift a{0) increases while the dilaton ip, the gauge field Q and 
the Higgs field / condense in regions closer to the hypermonopole core. 

In the equations of motion, n also affects the stress energy tensor produced by 
the Higgs and gauge fields. As earlier mentioned, the quantity S± in Eq. ( 134|) has a 
dependence in {n — l){n — 4) which implies the same negative contribution in seven 
{n = 2) and eight-dimensions {n = 3). However, the kinetic terms are multiplied by 
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Figure 3. Hypermonopole-forming field profiles in seven (n = 2) to ten {n = 5) 
dimensions. The coupling constants have been fixed to fiducial values a = 0.1, e = 0.2 
and /3 = 1. Gravitational effects are enhanced by increasing n with the exception of 
the metric coefficient uj(r). As can be checked in the lower panel, and in Eq. (j34p . 
deviations from fiat space are first reduced by increasing n from 2 to 3, and then 
amplified for n > 3. 

a factor {n — 2) and vanish in seven dimensions. The overall orthoradial energy has 
therefore a non-trivial dependence in the number of dimension for n < 3, but should 
become monotonic for n > 3. We can verify in the lower panel of Fig. [3] that deviations 
from fiat space in the metric coefficient u{r) are first reduced when n goes from 2 to 3 
before increasing again for n > 3. 

In conclusion, we have found hypermonopole solutions for all tested values of n. 
All the solutions exhibit the typical condensation of Higgs and gauge fields encountered 
in topological defect configurations. Since the dilaton also condenses, we do have an 
effective varying Planck mass between the brane and the bulk. There is also a strong 
gravitational redshift traced by the extra-dimensional profile of cr(r). Finally, at the 
matter field condensation radius, the metric coefficient u{r) remains almost constant and 
the extra-dimensions become cylindrically shaped. In the next section, we solve for the 
propagation of spin two fiuctuations inside such a background and show that gravitons 
becomes resonant on the length scales at which the extra-dimensions are strongly curved. 
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3. Resonant gravitons 

Restricting our attention to transverse and traceless four-dimensional tensor fluctuations 
h^^, the metric assumes the form given in Eq. f lT^ with the replacement rj^^, — )■ rj^y + h^y. 
Their linearised equation of motion can be obtained by expanding the action in Eq. ([3]) 
at second order and has already been derived in Ref . [31] . Defining the conformal radial 
distance z{r) and the rescaled tensor fluctuations by 

z = m^ f e-'^(")/MM, = e^'^/^+^/^zu'^/^h^,, (39) 

Jo 

the equation of motion for the gravitons reads [31] 



- C + + ^' - - a,j s,,y = 0, (40) 

where a prime stands for derivative with respect to z. The quantity W{z) is a 
superpotential given by 

W{z)^\a'^'^^ + U\ (41) 
4 Z w Z 



while A„ is the Laplace-Beltrami operator on the n-sphere S"" 

1 / 9^ n-i d 

1=1 ' 

and Qa is the four-dimensional d'Alembertian 



— 1/9^ n — id\ 



□4 = -\r]^iudf,dy . (43) 



In order to solve Eq. ( I40l) . it is convenient to perform a four-dimensional Fourier 
transform and a decomposition over the hyperspherical harmonics f^ such that the 
mode functions UM,e{z) satisfy 



Um p + 



um/ = 0. (44) 



W + W' + e{i + n-l)^ - M' 

w 

Here and — £(£ + n — 1) are the respective eigenvalues of the d'Alembertian and 
Laplace-Beltrami operator (in Higgs mass unit). This equation assumes the form of a 
Schrodinger equation of a supersymmetric quantum mechanical system [46j . The central 
potentials associated with fermionic- and bosonic-like excitations are given by 

V2{z) = W^ + W\ Vi{z) = W^-W'. (45) 

Omitting the tensor indices, the ground state .^0 of Eq- (1441) is the solution obtained for 
= i = and satisfies 

j.^w){-±^w)^,^0, (46) 
I.e. 

6 oc tu"/2e3-/4+^/2. (47) 
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This zero mode is not normalisable asymptotically. The ground state of the superpartner 
potential Vi is similarly obtained by swapping both terms in Eq. f H6|) and is given by 
^0 oc 1/^0- This time, it is not normalisable in the hypermonopole core for z ^ 0. As 
a result, there are no massless gravitons trapped on the brane and "supersymmetry" is 
broken by the solutions we are interested in (the spectrum associated with V2 and Vi 
do not match). Notice that the supersymmetric properties of Eq. (jUj) ensures that the 
spectrum is positive and no tachyonic propagation modes can be present. 

In order to solve Eq. ( I40l) in general, we assume that the mode functions um/ are 
normalised such that 



^/,^(^i)^M,^(^2)dM = 6{zi - Z2). (48) 
It is now straightforward to check that the Green function for ^^j, reads 







G^(Xi;X2) = - J 



d> 



3«Pm(^1-^2) 



X E fr^e^)frie.) [ .yf^'^^T' . m 

^ J M'' + p- [p" + le) 

Capital letters have been used for A^-dimensional coordinates, bold characters for the 
n-dimensional vectors lying on S^, and arrows for the usual three-dimensional vectors 
on the brane. Let us mention that we will not need to specify an explicit expression for 
the /J" functions and solely assume they form an orthonormal basis such that 

Y^frmirm = [jn(^i)]-^5"(^i - ^2), (50) 

where J'n{G)d9i ■ ■ ■ d9n is the infinitesimal surface element on the unit n-sphere so that 

X = sin'^-^(^i) sin"-2(02) . . . sin(ft„_i). (51) 

From the Green function, we can derive h^i, for any additional stress energy tensor 
on the brane. Considering an additional transverse and traceless four-dimensional source 
s^v{x) inducing a A^-dimensional linear stress-tensor perturbation of the form (in Higgs 
mass unit) 

5T^,{X) = \im^-^[Me)]-H^{0)s^,{x), (52) 

^->-o z 

the tensor fluctuations at X = Xi are given by Eq. ([8]) and read 

h (X.^- -3a(^i)/4 -V>(^i)/2 

W ' (Zij 

X y d^X2G'5(Xi;X2)tu"/2(^2)e''^(^^)/'e-^(^^)/25T^,(X2), (53) 

with = VTL^^'^H^ . The only unknowns are the mode functions UM/i^) entering the 
definition of the Green function in Eq. ( I49p . and solution of Eq. flH|) . It is instructive 
to solve them assuming no-dilaton and fiat space-time, i.e. without the presence of the 
hyp ermonop ole . 
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3.1. Flat spacetime 

Assuming ip = a = as well a.s w = z along the extra-dimensions, Eq. ( I44l) is a Bessel 
equations whose regular solutions in the origin read [17] 

uI,^, = Vm^Ju{Mz), (54) 

with 

(55) 



u = 

The orthonormalisation properties of the Bessel functions |48] automatically ensures 
that Eq. ( 1451) is satisfied. Plugging Eq. ( 154|) into Eq. ( 15^ and looking for solutions 
sourced by perturbations of the form fl52|) gives 



X 



(56) 



The function Dm{x) is the four-dimensional retarded propagator defined by 



d^p 



{2nY + {p)' - {p' 
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M 



2s 



(57) 



with 



(X 



0\2 



— (x)^. The term in Z2 shows that only the hyperspherical harmonics 



with zero eigenvalues i = contribute to the interactions sourced on the brane {z2 = 0). 
For static sources, performing the previous integrations and evaluating the solution also 
on the brane {zi = 0) yields 



2k' 



1/0 (0)1 



h 2"-i 



n + 1 



d^X2S.Jx2) X / dMM 



-M\Ax\ 

\Ax\ 



(5J 



where Ax = xi — X2. The last term in the previous equation is the Laplace transform 
of M" which is T{n + l)/|Af|"+^ From Eq. (ED]), one has 

-1 



|/o(0) 



Jn(^)d"^ 



1 

5" ' 



where 5" is the hypersurface of the unit n-sphere: 

2^(«+i)/2 



5" 



n + 1 



(59) 



(60) 



After having restored the dimensions, Eq. (!58|) simplifies into 



2k' 



{n + 2)S' 



s^,u{x2) 



\Xi - X2\ 



\n+2 



(61) 
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which is the standard hnearised solution of the Einstein equations in = 5+n spacetime 
dimensions. The iV-dimensional Newton constant also matches with the standard value 

= {n + 3)5"+^ • ^^^^ 
3.2. Inside the hypermonopole 

Inside the hypermonopole, the tensor fluctuations can be derived in a similar way. One 
should flrst keep the factors involving vj, a and if). In fact, as can be seen from Eq. ( 153|1 . 
by taking both the source and the observer on the brane, all factors involving a and 
w cancel, solely the dilaton rescales the gravitational coupling constant by exp[— '?/'(0)]. 
The mode function UM,i are no longer the same but for both the source and the observer 
on the brane, only their value in 2; = enters the calculation. Furthermore, since the 
only hyperspherical harmonic which is non-zero at = is /q, only the / = modes 
contribute to the tensors fluctuations. In fact, by defining the spectral density 

P(^) = — — 77^ ' (63) 



\u 



M 



one arrives at 

(n + 2)5"+^r(n + l) 
where we have defined the Laplace transform 

poo 

C^{q{M)}= / e-^^^g(M)dM. (65) 



9k-^ p-V'(O) r \ 



^0 

A four-dimensional behaviour can be recovered if the Laplace transform has a weak 
dependence in |Aa;|, which is precisely the case when gravitons become resonant with a 
mass m„. Taking toy example 

p{M) = QQ + Gm^5{M -m^), (66) 

where ^0 and G are two constants, one gets 

£|A.i{M>(M)} = ^j^^^ + Gml+^e-^^\'''\ (67) 



< mg|Af| < 1 (6^ 



whose second term dominates and is almost constant in the range 

- -p/ , ^\-|l/(n+l) 

G 

provided qq is sufficiently small. Outside of this range, the inverse power term dominates 
and the tensor fiuctuations are that of A^- dimensional gravity. The upper bound in Eq. 
fl68|) is satisfied for graviton resonances which are light enough, i.e. — > 0, whereas 
the lower bound requires a strong peaked resonance, i.e. a long lived graviton having 
C 3> f?o- In the following, we show that such a situation generically occurs inside the 
hypermonopoles: gravitons become strongly resonant due to the positive curvature of 
the extra-dimensions. 
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Figure 4. Superpartner potentials Vi and V2 for the eight-dimensional hypermonopole 
of Sec. [2] Although they do not support bound states, nietastable massive gravitons 
are trapped in the confining nest around z ~ 200. The peaks are associated with the 
regions of maximum curvature of Fig. [T] (same parameters). 



3.3. Graviton spectral density 

Since the extra-dimensions are asymptotically flat, the mode functions far from the 
hypermonopole core are of the form given by Eq. (15^ . Using these Bessel wave- form 
as asymptotic initial conditions, we have numerically solved Eq. ( I44ll in the background 
fields of Fig. [T] In Fig. HJ we have represented the superpartner potential Vi and V2 as a 
function of the conformal radial distance z. Solving for the mode equation in V2{z) gives 
the spectral density plotted in Fig. |3 For this configuration, we have observed more 
than fifteen trapped gravitons, the hghtest having a spectral density profile typical of a 
Breit-Wigner distribution 

PiM) ^ —j- ^ , (69) 

whose best fit gives 



7g 



mg = 2.4103240725 x lO^^mh, 
7g = 1.12 X lO^^mh 

I = 1.19 X 10^ (70) 

In order to properly resolve the width of these resonance, we have implemented a 
recursive local adaptive mesh refinement coupled to the more usual Runge-Kutta 
integration of Eq. (jH]). For such very thin resonances, the toy model of Eq. fl66|) is 
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0.001 0.01 0.1 



M/m, 

h 

Figure 5. Spectral density with respect to the graviton mass for the eight- 
dimensional hypermonopole. The resonances correspond to the massive metastable 
states associated with the potential V2{z) plotted in Fig. ID More than fifteen gravitons 
end up being trapped. The constant behaviour at small and large distances shows that 
gravity is iV-dimensional in these regions, albeit with a different effective Newton 
constant due to the gravitational redshift. 



a good approximation for M ~ mg, the constant C being given by the integral of 
Eq. ([69]). One finds 

C = (n + 2 arctan ^) ~ tt/^ , (71) 

2mg V IgJmg 

where the last expression is accurate only for 7g ^ nig. For the best fit values of Eq. (I70|) . 
one finds C ^ 1.7. As can be seen in Fig. (|5]), for light masses p{M) = Qq = 8.2 x 10~^ 
such that the lower bound in Eq. f pS]) is about 10^^. We therefore expect this resonance 
to change the standard eight- dimensional gravity law on distances covering not more that 
two-orders of magnitude around the scale l/rUg, which is far to short to be interesting 
for cosmological purpose. 

However, we do not see any reasons preventing the existence of cosmologically 
interesting solutions, i.e. much lighter gravitons. Indeed, increasing a appears to push 
up the potential barrier in Fig. H] whereas reducing e increases the width of the potential 
barrier. Small values of e have the effect of delocalising the gauge field and this ends up 
spreading its energy density over the extra-dimensions. Both of these parameters could 
therefore be somehow adjusted to obtained much lighter graviton resonances. As the 
numbers reported in Eq. (170|) suggest, the precise determination of lighter resonances is 
made difficult due to numerical limitations, the machine precision accuracy not covering 
more than 16 orders of magnitude on usual computers is already saturated by 7g// in 
Eq. dZOD. 
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10° lo' 10' 10^' 10"' 10° lo' 10^ 10^ lO"* 10^ 



X m, X m, 

h h 

Figure 6. Rescaled Laplace transform in linear (left) and logarithmic scale (right) 
associated with the spectral density of Fig. [5] and plotted as a function of the distance 
to the source. The Newton law is eight-dimensional when this quantity is constant as it 
occurs on small and large distances. The lightest graviton resonance of Eq. ([70]) reduces 
locally the dimensionality of the Newton force by a fractional power of exponent 2.4 
(best fit). The strong variation in amplitude from small to large scales comes from the 
hypernionopole induced gravitational redshift. 



3.4- Deviations from Newton 

3.4-1 ■ Dimensional reduction In order to complete the discussion of the previous 
section, we have computed the Laplace transform directly from the spectral density 
found in Fig. |5l Fig. |6]shows the rescaled quantity \x\^^^ C\^\{M'^ p{M)} as a function of 
the distance to the source. Gravity is A^- dimensional when this quantity is constant as it 
occurs at small and large distances. The strong variation in amplitude, also visible in the 
smooth change of the spectral density in Fig. [5l comes again from the hypermonopole 
gravitational redshift (see the a profile in Fig. [1]). The graviton resonance of Eq. (ITOll 
is responsible of the peak located around x ~ l/'^g- A best power fit of the potential 
nearby this region shows that the Newton law is dimensionally reduced to This 
is not yet a four-dimensional Newton law due to the previously discussed numerical 
limitations to obtain a light enough resonance. 



3.4-2. Effective gravitational coupling From Eqs. (!62l) . (!64l) and (!67ll . one can extract 
the effective Newton constant G which would be measured by a four- dimensional 
observer in the three expected regimes. At very small distances. 
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the spectral density is constant p(M) = 1 and thus gravity is seven dimensionaL 
The measured Newton constant is however reduced by the dilaton condensation in the 
hypermonopole core: 

Go = e-^(°)G^ . (73) 

In the intermediate regions, those verifying Eq. f l68l) . gravity is driven by the metastable 

hghtest graviton and the Newton law is four- dimensional with an effective Newton 

constant given by 

^ ^^(o)giv n + 3 G 

2 n + 2r(n + l) s ^ > 

This equation makes clear that a light graviton, required for the dimensional reduction 
of the Newton law, will necessarily induce a small effective Newton constant thereby 
addressing the mass hierarchy problem. In terms of the reduced Planck masses, Eq. fTMl) 
can be recast into 

,n + 2r(n + 1) M"+3 
n + 3 C ^ 



Mj, = 2e««)!i±iiA!i±ii^, (75) 



where M|i = I/G4 and M^-^ = I/Gn- 
Finally, on the largest length scales, 

|Af|mg > 1, (76) 

gravity becomes again A^- dimensional but with a much weaker Newton constant since 
now p(M — )■ 0) = ^ 1- The measured Newton constant is now given by 

Goo = f?oe-^(°)G^ = fPoGo « G,. (77) 

As a numerical application, we can determine the order of magnitude of the A^- 
dimensional Planck mass such that the graviton mass is of the same order than the 
cosmological constant energy scale, i.e. ~ 10~^^GeV. From Eq. (175|) . one gets 

1 l/(n+3) 



M 



T(n + 1) 



25-lln 



X lO^^GeV, (78) 



which is down to the 10 GeV scale already for n = 2. Notice that the lowest scale 
at which A^- dimensional gravity shows up is fixed by the value of Qq and not by nig. 
This quantity coming only from the gravitational redshift ct(0), it can actually be made 
arbitrarily small for order one coupling constants, i.e. for Higgs vacuum expectation 
values also around the 10 GeV scale. In fact, as suggested by the previous equation, the 
mass hierarchy mechanism advocated here is so efficient that it has a natural preference 
for very small iV-dimensional Planck masses. 

4. Conclusion 

In this paper we have shown that metastable massive gravitons generically exist in the 
four- dimensional core of any self-gravitating hypermonopoles formed by the breakdown 
of an S'O(ne) symmetry in (rZe + 4) dimensions, provided Ue > 3. 
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Since the extra-dimensional spacetime is of infinite volume and asymptotically fiat, 
these resonances induce a DGP-like gravity confinement mechanism in the core. For 
light enough resonances, gravity is A^-dimensional at small and large distances, but can 
be four-dimensional on some intermediate range. The numerical determination of such 
a light and long-lived resonance may be however a non-trivial problem due to finite 
numerical accuracy. Moreover, we have shown that, in this regime, the effective four- 
dimensional Planck mass is proportional to an inverse power of the graviton mass; this 
one being extremely light, the mass hierarchy problem ends up being naturally addressed 
in our setup. The strong decay of the gravity law at large distances, coming from both 
the higher-dimensionality and the strong gravitational redshift, might be of interest to 
explain the current cosmic acceleration. 

Finally, we would like to emphasize that these models still remain unexplored on 
various aspects which may compromise, or not, their viability. Here, we have only solved 
the propagation of spin two fluctuations which decouple from the background fields. The 
model has however vector and scalar modes which may propagate and might also be 
confined in the core. Solving for their propagation is a challenging problem since they 
will be necessarily coupled to all of hypermonopole-forming vector and scalar fields. We 
leave the second order perturbation of Eq. ([3]) in the scalar and vector modes for a future 
work. 
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5. Appendix 



In this appendix, we present the Einstein tensor Gab and the stress tensor Tab- Inserted 
into Eq. (IHD, we then arrive at the equations of motion for the gravity sector ( 126|) . ( 1271) 
and (l28|l . For the metric given in Eq. (IT^ . the Einstein tensor is given by 



G 



-a + -a^ +n h 

2 2 uj 



zu n{n — 1) f w\'^ 3n&w n{n — 1) 



VJ 



2 U7 



2^2 



Gf 



ml 



2^ ^ 2 



vj\ „ avj nin — 1) 
- +2n 



Geidi - mlge^e. 



+ 2(n-l] 



5 , , o zu 
I w 

aw {n — l){n — 2) f roX^ {n — l){n — 2) 



zu 



zu 



2vu^ 



(79) 



With the Higgs and gauge fields specified by Eqs. ffT7|) . and fl?T]) . the stress 
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tensor (ITU]) becomes 



T, 



T 

J. r", 



Tf>, 



2 2 

V m 



-f-n 



fQ 



n -2 n{n — Ij 
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